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I. INTRODUCTION 

Among the different approximate analytical ap- 
proaches to the calculation of the Casimir energy and 
related objects [T], a particularly interesting one has been 
developed some time ago [2], for cases where the mirrors' 
configurations can be described as a small geometrical de- 
formation of the simplest possible case: two infinite, fiat, 
parallel mirrors. We shall focus here, in particular, on 
cases where just one of the surfaces defined by the (zero 
width) mirrors is deformed. Since one of the surfaces is, 
indeed, flat, it is always possible to adopt Cartesian co- 
ordinates xi,X2,X3 in such a way that surface coincides 
with the = plane. The other mirror's surface will, 
in turn, be assumed to be a small deformation of the 
X3 = a (a > 0) plane, namely: X3 = a + r)(xu). Here, 
r\ is a scalar function of X11 = {x\,X2), coordinates which 
are then parallel to the flat mirror surface. Hence, 77 
represents the (assumed) small departure from a config- 
uration corresponding to two fiat parallel mirrors (i.e., 
x 3 = 0,a). 

As already mentioned, for this, and similar situations, 
a perturbative expansion in powers of //(xy) has been de- 
veloped, in order to compute the Casimir energy, under 
the assumption that |j7(xm)| << a. This program has 
been carried out up to the second order in r\ [2], and 
applied, for example, to the evaluation of the effect of 
corrugations on the Casimir energy [5H3] . The perturba- 
tive calculations of Refs. [5J [1] are based on a functional 
approach, that we will also follow here. On the other 
hand, in Ref. [3] a scattering approach is used, and a 
general formula for the Casimir energy is given in terms 
of the reflection coefficients associated to the mirrors. In 
this approach, one can incorporate finite temperature, fi- 
nite conductivity and roughness corrections. However, as 
far as we know, all previous calculations for the normal 
force [2-4; have been performed up to second order in r\. 

In recent years, accurate numerical evaluations of the 



Casimir energy for several geometries, including corru- 
gated surfaces of arbitrary amplitude, have become avail- 
able. In spite of that, approximate analytical results are 
always welcome, for a variety of reasons. Among them, 
to improve the physical understanding of the dependence 
of the Casimir force with the geometry, as well as to use 
them as benchmarks of complex numerical calculations. 



It is the aim of this paper to present new analytical 
results for the perturbative expansion, which had been 
implemented in previous works up to the second order 
in rj, by calculating the two subsequent corrections, of 
order three and fourth in the same expansion (note that 
a similar program has been carried out in Ref. [5], to 
compute the lateral force between corrugated surfaces up 
to the fourth order in the amplitude). 



The structure of this paper is as follows: in Sec. |TIJ we 
introduce the system considered and derive the (formal) 
expression for the Casimir energy, in terms of an "effec- 
tive action" , using a functional approach. In Sec. |III[ we 
introduce the perturbative expansion for the effective ac- 
tion, focusing on the expressions corresponding to terms 
of up to the fourth order. In Sec. \YV\ we present the ex- 
plicit results up to the fourth order term. In Sec. [V] we 
analyze further approximations and resummations of the 
perturbative expansion. In particular, we will see that 
the proximity force approximation (PFA) can be under- 
stood as a resummation of the perturbative expansion 
when the form factors involved in the expansion are re- 



placed by their zero-momentum values. In Sec. VI we ap- 
ply the results to compute the correction to the Casimir 
energy up to the fourth order for a sinusoidally corru- 
gated surface, and up to the third order for a corrugation 
described by a combination of two sinusoidal functions 
of different wavelengths. We summarize our findings in 
Sec. EED 
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II. THE SYSTEM 

We adopt Euclidean conventions, whereby spacetimc 
coordinates are denoted by x^ = x^ (ju = 0, 1, 2, 3), 
x being the imaginary time and Xi, (i= 1,2,3) the 
spatial coordinates. We shall also use the notation 
x\\ = (xo, xi, X2) = (xo, xy), which shall be used as space- 
time coordinates on each mirror. 

As mentioned above, the scalar field p satisfies Dirich- 
let boundary conditions on two surfaces: L, which is the 
x 3 = plane, while the other, R, can be represented by 
a single function of xm : 



L) x 3 = 

R) x 3 = a + ?7(x||) 



(1) 



Following the functional approach to the Casimir ef- 
fect, we introduce iJ, the zero temperature limit of the 
partition function for ip in the presence of the two mir- 
rors: 



T>tp 5 L (cp)S R (cp) e~ So ^ , 



(2) 



where So is the massless real scalar field Euclidean action 



s (v) = - I d 4 x(d<p) 2 



(3) 



while 5l (Sr) imposes Dirichlet boundary conditions on 
the L (R) surface. Because of several reasons, it becomes 
convenient to generalize, at this point, the boundary con- 
dition on the R mirror, by assuming that 77 may also de- 
pend on Xq, the Euclidean time. Namely, ?y(x||) — > ?7(a;||). 
One of the reasons for doing this is that the results thus 
obtained become more symmetrical, in the sense that the 
coordinates afn shall appear on an equal footing. Besides, 
as a byproduct, one could make use of those results in 
dynamical Casimir effect situations. Note that, strictly 
speaking, Z will not be a quantum partition function 
when 77 depends on time. We do keep, however, the same 
notation (Z) for the resulting object. Static Casimir ef- 
fect results will be obtained by simply dropping the time 
dependence of rj, at the end of the calculation. 

Finally, note that Z, evaluated with a time dependent 
r], may also be interpreted as the classical partition func- 
tion for a static system in 4 spatial dimensions, in the 
presence of two (static) boundaries. 

To proceed, note that the two functional delta func- 
tions may be exponentiated by introducing two auxiliary 
fields, Al and A^, to obtain for Z the equivalent expres- 
sion: 



Z = VpVX L VX R e 



-S(<p;\l,>>r) 



(4) 



with 



S((p;X L ,X R ) = S (<p) -i J d 4 x[X L (x\\)S(x 3 ) 
+X R (x l] )8(x 3 -<>p{x l] ))]<p{x) , (5) 



where ijj(x\\ ) = a + r/(x\\). 

Integrating out tp, we see that a Zq factor, the parti- 
tion function for the field tp in the absence of boundary 
conditions, factors out, while the rest is an integral over 
the auxiliary fields: 



Z = Zn 



VX L VX R 



g-2 S d3x \\ I d3 V\\ E«,# ^c(x ]] )T afl (x ]] ,y ]] )X l3 {y^) ^ ^ 

where a, j3 = L, R and we have introduced the object: 

T LL (x {] ,y {] ) = (x^OK-d 2 )- 1 ^) (7) 

T LR (x hyil ) = (x^^K-d 2 )- 1 ^,^)) (8) 

T^ix^) = (x\ { ,fP{x\ { )\(-Sf , )- 1 \y\ { ,0) (9) 

T RR (x {] ,y {] ) = (ar||,^(*||)|(-c> 2 )- 1 |</|| ) V>a/||) . (10) 

We have used a 'bra-kef notation as a convenient way 
to denote matrix elements of operators, d 2 is the four- 
dimensional Laplacian. Thus, 



(zK-a 2 )- 1 ^) = 



d 4 k e ik < x -^ 



(27T) 



k 2 



(11) 



The four kernels T a( g can be naturally regarded as matrix 
elements of a 2 x 2 kernel matrix: T = (T a p) a /3 _ L R - 

The vacuum energy of the system, £v ac , obtained by 
subtracting the zero-point energy of a free field with triv- 
ial boundary conditions at infinity, may then be obtained 
as follows: 



lim 

T->oo 



? /=»)( x ii) 



(12) 



where T is the extent of the time dimension, and T is 
given by: 



r = -io g (- 



= ^Tr(logT) 



(13) 



wich resembles an effective action in the presence of a 
background field r](x\\). Note that the trace above is 
meant to act on both discrete and continuous indices. 

The Casimir energy is then obtained from the vacuum 
energy E vac , by simply discarding contributions due to 
the vacuum field distortion produced when the mirrors 
arc infinitely far apart. This 'self-energy' piece, is irrele- 
vant for the kind of physical observable we have in mind 
here, and shall therefore be subtracted. 

In the next section, we construct the expansion of T in 
powers of 77, whence the information about the Casimir 
energy will be extracted from (12), after discarding self- 
energy terms. 



III. PERTURB ATIVE EXPANSION 

The procedure to obtain the formal expression for the 
expansion of T is rather straightforward. Indeed, assum- 
ing that denotes the term of order i in the expansion 
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of T in powers of 77, 

T = T (0) + T (1) + T (2) + T (3) + T (4) + . . . , (14) 
we obtain an expansion for F of the form 

r = r (0) + r (D + r (2) + r (3) + r (4) + _ (15) 

Terms of up to the second order for T are then given by 
r (°) = lTrlog(T(°)) , r« = ^Tr [(Tt ))- 1 !^] (16) 



and 



where: 



p(2) = p(2,l) + p(2,2) 



(17) 



p(2,l) = i Tr 

2 

p( 2 < 2 ) = _Itt 
4 



/ T (0)\-1 T (2)' 

(■f(0))-l T (l)( T (0))-l T (l)l _ ( 18 ) 



Collecting all the terms of third order in 77, we see that 
I^ 3 ) may be expressed as the sum of three contributions: 



p(3) = p(3,l) + p(3,2) + p(3,3) 



where 



r (3,i) = 4 Tr 
2 



^ T (0)^-l T (3) 



(19) 
(20) 



p(3,2) = _^ Tr J( T (0))-l T (l)( T (0))-l T (2)j ; (21) 



and 



r (3,3) = ^ Tr ^ T (0)^-i T (i)^ T (0)^-i T (i)^ T (0)^-i T (i; 

6 '(22) 
Finally, r' 4 ) is the sum of five independent terms, 
which in turn result from collecting all the terms of fourth 
order in the expansion for T: 



p(4) = ^r (4 < 



i) 



(23) 



and 

p(4,5) 



(28) 



Tr 



^(0)^-l-jp(l)^(0)-j-lj(l)^(0)-j-lqp(l)^(0)^-l-jp(l) 



To proceed to the evaluation of the explicit form for 
each term, it is convenient to introduce ?7(fe||), the Fourier 
transform of T](xu), 



xn e 



ik 1 -x 



V{x\\) ■ 



(29) 



Also, we note that (TH -*) 1 , which is translation invari- 
ant, may be written as follows: 

(T(°>rV||,Z||) = (T (0) )-V||-Z||) 

' rkn e^r^u-xdom) , (30) 



(2tt) 3 



with: 



_ e -|fe|ik 



^ = T^^(. e \la ) • (3D 

Another important ingredient is the form of the matrix 
elements of T to the desired order. Regarding the diag- 
onal elements, one sees that T ( ll = for j > 0, while 

(i) 

Trr = for odd j. Thus, for the diagonal elements, we 
need: 

(32) 

and 

(33) 

For the off-diagonal matrix elements, their general form 
at an arbitrary order may be written in a rather compact 
form. Indeed, 



with: 



p(4,l) = ^ 

2 



( T (0))-1 T (4) 



(24) 



(27T) 



(34) 



p(4,2) = _ 4 Tr 

2 



( T (0))-i T (i)( T (0))-i T (3)l ) (- 25 ) 



(35) 



p(4,3) = _£ Tr 

4 



(T (0) r l T (2) (T (0) r l T (2)j ) (2g) 



IV. RESULTS 



p(4,4) = ^r^ 
2 



( T (0))-1 T (1)( T (0))-1 T (1)( T (0))-1 T (2) 



(27) 



Although the terms of order zero, one and two arc 
already known, we present, for the sake of completeness, 
the results corresponding to those orders as well. 
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A. Zeroth order 

The zeroth-order term is simply obtained by recalling 
the form of T^°\ from which one just needs to subtract 
its a — > oo limit. After that subtraction, one obtains a 
result that may be written in the following way: 



lim 

T->oo 



r(°) 



T 



l2 

2 



d 3 P\\ 



log(l 



-2p\\a\ 



(36) 



where L is a length that measures the size of the mirror 
along the parallel directions. 



Thus, the Casimir energy due to this term, E^ c , 
comes: 



be- 



L 2 



d 3 p\\ 



log 1 



-2|p|||o"V _ 



TT 2 L 2 

1440 a 3 



as expected. 



(37) 



First order 



It is quite straightforward to compute the first order 
term I^ 1 '. Introducing a function B, defined as: 



B(q\\) 



o 2 l9|||o 



(38) 



where q» is a 3-vector, we obtain 



d 3 pn 



r (1) = / S fl (Pii)hl x / ^ii^ii) 



d 3 xu r](xu ) 



480a 4 



(39) 



Since, by definition, r\ measures the departure from a 
flat parallel mirrors configuration, one could impose on it 
the condition that the integral / d 3 x\\ r\{x\\) equals zero. 
Were it not the case, this condition could nevertheless 
have been achieved by subtracting a constant (which has 
to be added to a) from r\. Thus, in this case rW = 0. 



C. Second order 

The second order result has been known for quite some 
time. In Fourier space, it may be written as follows: 



r (2) = 



with: 

/ (2) (fc 



1 

2 

d 3 



d 3 k n 



kf {2) (h)\v(k\\)\ 



(27T) 



(40) 



[B( P ||) J B(p||+fc|,)+ B(p l]+ 

(41) 



(2tt) 3 

X |p|| +fc||| 

From the structure of p0| ), it seems that, for the pertur- 
bative expansion to make sense, an extra necessary con- 
dition (besides \r]\ << a) is that the integral over fen has 



to be well-defined. However, we can, and will, consider 
cases where rj(k\\) is a generalized function. In particu- 
lar, <5- like, concentrated around a particular value of ku ■ 
In this case, all the correction terms (not just the second 
order one) become proportional to the size of the system, 
due to the periodicity of the perturbation. Of course, the 
same is always true of the zeroth order term. Thus, even 
for this singular case, one could make sense of the pertur- 
bative corrections, at least if there is a region where the 
corrections are reasonably small, in comparison with the 
zeroth order term, after factorizing the (common) spatial 
size factor. 

On the other hand, when 77 (fey ) is a regular function, 
a nontrivial condition emerging from ( 40 ) proceeds from 



requiring its convergence at large momentum. It may be 
seen that, at second order, this amounts to: 



d 3 fc|| 



\h\Uh)\ 



< oo 



(42) 



In a realistic situation, the Fourier spectrum of the de- 
formation rj should have a cutoff; for instance, a contin- 
uous description of the mirrors should, at some point, be 
replaced by a discrete (lattice) one, which introduces a 
large momentum cutoff of the order of the lattice spacing. 



D. Third order 

The contribution denoted by r*- 3,1 ) is ultralocal, under- 
standing by that an integral of a term involving (three) 
?7's at a single point (without derivatives). The first ex- 
ample of an ultralocal contribution was the first order 
term, but it could, as we have seen, be assumed to be 
equal to zero. On the contrary, r' 3,1 ', given explicitly 
by: 



p(3,l) = 



1 



d 3 P|| 
3! J J2nf' 



B(p||)b||| 3 x / d 3 x {] [r,(x l{ )] 



1 3 _ T-( 3 ) 



(43) 

(3) 

does not necessarily vanish. We have used T\ , to denote 
this contribution, the only ultralocal term at the third 
order. 

On the other hand, r^ 3,2 ) is 'bilocal', i.e., it is cubic 
in rj and involves products of r\ at two (eventually equal) 
points. It is the only bilocal piece of T^ 3 ). Its explicit 
form is 



r( 3 - 2 ) 



d 3 fc,i 



(27T) 



/f ) (fc||»|lh 2 (-/qi) 



3 J b 



.(3) 
b ' 



( 44 ) 

(note the Fourier transform of r\ ), where we introduced 

(3) 

j h , the bilocal kernel at the third order: 



(2tt) 3 



b (p\\) hilhi + h 



(45) 



Finally, r^ 3 ' 3 ) is trilocal, since the three 77's appear at 
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(generally) different spacetime points: 



r( 3 . 3 ) 



f 



3! J (2tt) 3 (2tt) 



x rj(-k {] - q {] ) = if > , 
where the /( 3 ) kernel is given by: 



(46) 



/* (3) (Vll) = 2 



d 3 pii 



|*|| + i|| +P\\\ + P\\ 



(2tt) 3 

x 4B(fc|| +l\\ +P\\)B(l l{ +p {l )B( Pll ) 

+ B(k\\ +P||)-B(i|| 

+ 5B(fc|| + /y + V\\ ) B (P\\ ) + 3B(fc|| + Z 



(47) 



The full third order term is then given by the sum of 
the previously derived local, bilocal and trilocal contri- 
butions: 



r (3) 



,(3) 



r (3) , r (3) 
1 b "T 1 t 



(48) 



E. Fourth order 

As we have just done for the third order contribution, 
we present the partial results corresponding to what we 
denoted as r^ 4 '^, for the different values of Z, which run 
from 1 to 5. 



The I^ 4 ' 1 ) term, given by (24), contains both local and 



bilocal parts, the local one being a-independent. After 
discarding the a-independent part, the explicit form of 
the term becomes: 



r (4,l) = 

with 

/ (4 ' 1} (A||) 



d 3 k 



(2tt 

V 3 (-k\\)v(k\\)] 



d 3 P\\ 
(2nY 



B(p\\)\p\\\\p\\ 



(49) 



(50) 



The next term, I^ 4,2 ) may be written as follows: 
1 r d 3 h 



p(4,2) 



(2tt) : 



-J {i,2) (h)v 3 (-h)v(k\\) ■ (51) 



with 



/ (4,2) (*||) = -|/ (4,1) (fc||) 



4 r 

2 f d 3 pn 



(52) 



B(p» +k«)B(p«) \p»\\p» +k« 



3 J (2tt) 3 

We now write the result corresponding to T^ 4:3 \ which 



involves three different kernels: 



r( 4 - 3 ) = 1 
4! 



d 3 k\\ d 3 q\\ d% (43) 



(2tt) 3 (2tt) 3 (2tt) 3 

x vinWhW-h -q\\ - 

+ 3! / (2- (»■< 
x ?7 2 (— fcii — gn) 

1 r d 3 k 
+ 2! J (2tt) 



i yt' 3) (h)v 2 (h)v 2 (-h). (53) 



After some rather lengthy algebraic manipulations, we 
see that the form of the kernels introduced above is the 
following (again, discarding a-independent terms): 

d 3 p« 



f^ 3 \k h q h lii) = -6 



+B(fc|| +q\\ +/|| + p\\)+B(l\\ +py)] 
x |fc|, +q l{ +/|, + P| |||g|,+f|| +P|||b||l . ( 54 ) 

d 3 pii 



/ 2 (4 ' 3) (fcn,9 
and 



(2tt) 3 



/f 3) (fc||) 



B(p||)|/c||+g||+ P ||| 2 |a||+p|| 



(55) 



(56) 



Let us present now the results corresponding to the 
p( 4 > 4 ) contribution. We have found that it contains both 
4-local and 3-local terms: 



r (4,4) 



1 

4! 



d 3 k\\ d 3 q\\ d 3 l l{ 



ft A \h,q ll ,l ll )rj(k ll ) 



(2tt) 3 (2tt) 3 (2tt) 3 

x v(q\\)v( l \\)v(- k \\ - ?ii - 

x ?7 2 (-fc|| -q\{) , 
with the kernels: 



(57) 



-12 



d 3 p|| 

(2^y 



|A)|i +qu +L +p\\\ 



x +P||lb||| 

x 4B(*|| +g|| +P\\)B( qil +/|| 

+ &B(fc|| +/|| +P||) 

+ B(?n+Z||+p||)B(Z||+p||) (58) 

+ q\\ + l\\ +p\\)B(l\\ +p\\) + B(q\\ + l\\ +P||) , 

and 



# 4) (fc||,an) = 3 



d 3 p ii 



\k\\+q\\ +P\\\k\\ 



(2tt) 31 

W(k {{ +g,| +P\\)B{ qil +p,|) 
B(q« + pii)B(pii) + B(q« +pn 



(59) 
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The r< 4 ' 5 ) term is 4-local: 



p(4,5) = 



4! 



d 3 kn d 3 qn d 3 ln s 



x v(Q\\)v{ l \\)v(-h - q\\ - in) 



(60) 



with: 



/( 4 ' 5 )(fc||,q||,/||) = -6 



rf 3 



P\\ 



(2tt) e 



fc|| +q|| +/|| + ; 



x +p||||Z|| +p||||p||| 

x 8B(fc|| +g|| +Z|| +p||)B(g|| +/|| +P||)-B(/|| +P\\)B(p\\) 

+ 5B(fc|| +g ( | +/|| +P\\) 
+ 2B(k\\ +q\\ +/|| +P||)S(P||) 
+ B(fc|| + 9 || +/|, +p\\)B(ln +pu)B(p l{ ) 



+ B(k\\ +q\\ +/|| +P||)B(/|| +pn) 



(61) 



V. LOW AND HIGH MOMENTUM 
EXPANSIONS: CONNECTION WITH THE PFA 

In order to discuss the forthcoming approximations, we 
find it convenient to first write the previously considered 
perturbative expansion as follows: 



7T 2 L 2 

1440a 3 ! 7 



EL 



11 r d 3 k\p d 3 fc, W 
f ^ a 3+n J ~{2n) 



(2tt) 3 



x S{kf + ... + kf>)hM(akp, akf>) 



(62) 



The explicit form of the form factors l~S n \ up to n = 4, 
is completely determined by the results obtained in the 
previous section. Indeed, one just needs to include the 
proper S function factors in some of the terms we cal- 
culated, before adding them. In the equation above 
we have made explicit the fact that the (dimension- 
less) functions h^ n ' depend on the dimensionless variables 
(?) 

aK , 3 = 1) •••>«• 



Note that ( 62 1 has the structure of the general Taylor 



expansion of a functional in terms of its argument, an 
expression entirely analogous to the one we could use for 
the expansion of the effective action in a 2 + 1 dimen- 
sional quantum field theory, in terms of the 'field' rj, the 
form factors h^ n ' playing the role of the n-point vertex 
functions. 



A. PFA and the low momentum expansion 

Let us assume here that the function T?(xn ) is slowly 
varying. In terms of its Fourier transform, it will be 
peaked at zero momentum; therefore we can approximate 



M n )(a/4 1 \...,a^™ ) ) ^ /i (n) (0, ...,0). As a consequence: 



EL 



1440a 3 



n>l 



r ft< B >(0,..,0) / dWxn) 



In principle, one could evaluate the first terms by using 
the results of Sec|IV[ However, there is also a shortcut: 
for a constant 77 (xi 1 ) = 770, the vacuum energy is simply 



given by Eq.(63) with the replacement J ePxy ??(x||)" — s- 
L 2 r]g . But for this particular case we know, of course, 
that the answer is the Casimir energy between paral- 
lel plates separated by a distance a + 7/0 ■ Therefore, in 
the low momentum approximation the perturbative se- 
ries can be summed up, the result being 



1440 



d 2 xn 



+ ^(X||)) ; 



(64) 



which agrees with the PFA. 

Note also that the above suggests a nontrivial consis- 
tency check for our calculations: indeed, one can compute 
the first terms of the series in Eq.(63) using the explicit 



expressions given in Sec |IV| for the different contributions 
to the effective action r^ ,J K We present some details of 
that calculation in the Appendix. The result is 



1440a 3 



d 2 xi 



6( 



^2 



10( 



^3 



15 



a / 



(65) 



that agrees with the expansion in powers of rj of the PFA 
result. 

One could also go beyond the PFA, by expanding the 
form factors around fcii = 0. In this way, one should be 
able to recover the derivative expansion for the Casimir 
energy that we introduced in Ref. [6 . Once more, it is 
useful to have in mind the analogy with the usual expan- 
sions of the effective action in quantum field theory: the 
expansion in powers of the field corresponds to the per- 
turbative expansion in powers of 77, while the derivative 
expansion corresponds to an expansion in derivatives of 
tp = a + 77, i.e. the improved PFA [BJ. 



B. High momentum expansion 

We shall now consider the opposite limit, namely, 
|fc|||a ^> 1, in which the scale of variation of the shape 
of the surface is much shorter than the mean separation 
between surfaces. 

Let us begin by considering the second order term. In 
order to obtain t he la rge-momentum behavior of the form 
factor /( 2 ) in Eq.(41 ), we note that the first term is expo- 

-2o|fe||| 



nentially suppressed, i.e., it can be written as |fcii e 



7 



times a convergent integral. After a shift in the integra- 
tion variables, the second term can be approximated by 

d 3 p» 



f {2) (h) 



-2\h 



(2tt) 3 
d 3 p 



(2tt 



b(p\\) bnibu 
1 3 B (p\\)\p\\ 



fen 



240 a 4 



Inserting this result into Eq.(|40|) we obtain 

-2 



r (2) 



7T 

'480a 4 



\h\mh)\ 



(66) 



(67) 



Unlike the low momentum expansion, the result is a non- 
local functional of the shape of the surface. It can be 
rewritten in configuration space in terms of the nonlocal 
operator (—V 2 ) 1 / 2 as follows: 



r( 2 ) 



48^ /^ll^( x ll)(-V 2 ) 1/ ^(x|| 



(68) 



We now consider the third order terms. From Eq.(45) 
it is quite straightforward to see that 



if 3 W 



480 a 4 



(69) 



Note that this contribution grows quadratically with the 
momentum and becomes local in configuration space: 



r( 3 . 2 ) 

T 



7i~" 

960a 4 



d 2 xn ?7 2 (x||)V 2 ?7(xjj 



(70) 



The trilocal contribution is determined by the form factor 
given in Eq. ( 48 1 . For high momenta, the leading contri- 



bution comes from the last term, wich is linear in B. We 
obtain 



ft i3 \k h h) 



7T 2 |fc||(fc||+/||)| 



(71) 
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which produces a nonlocal contribution to the vacuum 
energy 



r( 3 > 3 ) 

T 



480a 4 
x (-V 2 ) 1/2 ??(x||) 



d 2 X|| 7 ? (x||)(-V 2 ) 1/2 f? (X||) 



(72) 



The fourth order terms can be treated in a similar way. 
One can show that, to leading order, they grow cubically 
with the momenta. 



VI. EXAMPLES 

In this section we present some applications of the per- 
turbative results obtained so far. As we will consider 
surfaces with periodic perturbations, the vacuum energy 
will be proportional to the size of the system. There- 
fore, we will compute the vacuum energy per unit area 




FIG. 1. Plot of the function gi defined in Eq.(73l, as a func- 
tion of the dimensionless variable x = q 1 a. Linear behavior 
of <72 for large x is shown. 



A. Sinusoidally corrugated surface 



We will first consider the simplest case of a sinusoidally 
corrugated surface ?/(x||) = esin^xi), with e <C a. For 
this particular corrugation, the first and third order cor- 
rections do vanish, so on general grounds we expect 



1440a 3 



1 +92 (qi a) 



(73) 

for some functions g 2 and (74. The function g 2 has already 
been calculated in Ref.[2]. For the sake of completeness, 
we present a plot of its numerical evaluation in Fig.l. In 
the zero momentum limit, the numerical results repro- 
duce the expected value, given in Eq.([65|. Moreover, the 
behavior for low momentum is quadratic [2j [6] , which is 
consistent with the existence of a derivative expansion for 
the Casimir energy [BJ. On the other hand, in the limit 
qia ^> 1, the plot reproduces the linear behavior g 2 ~ qia 



with the appropriate coefficient predicted in Eq.(67). 



The behavior of is shown in Fig. 2. Once more, in the 
zero momentum limit the numerical result is consistent 



with Eq.(65) (see inset in Fig. 2 for the small momentum 
behavior). On the other hand, in the high momentum 
limit, 34 ~ (gia) 3 . It is worth to note that the per- 
turbative expansion breaks down for very high momenta 
(short wavelength of the corrugations). Indeed, the ratio 
of the fourth to second order corrections is proportional 
to (egi) 2 . Therefore, the fourth order correction becomes 
more relevant when approaching the short wavelength 
limit. In Fig. 3 we plot the ratio between the fourth and 
second order correction coefficients, 34/(72: as a function 
of x = q\a. 




FIG. 2. Plot of the function 34 denned in Eq.([73f , as a func- 
tion of the dimensionless variable x = q±a. The inset shows 
the low momentum behavior of g^. 
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FIG. 3. 34/32 as a function of x. Fourth order correction be- 
comes more relevant when approaching the short wavelength 
limit. 



ily evaluated and gives 
2 

2 

1440 a 



^=Je(^) 2/(2)(9(1)) 

1 

^±{^) 2 9^ ( M (75) 



where — (0,q^,0, ..0). This is just the sum of the 
energies of each Fourier component of the correlation. 
In order to evaluate the third order contribution to the 



Casimir energy, recalling Eq.(48) we write 



p(3) _ f(3,l) , F (3,b) , o{3,t) 



(76) 



The different terms in the above equation can be evalu- 
ated computing the Fourier transform of rj(xu ) and using 



Eqs.(43l, (44|, and (461. Assuming that q[ 2) > q[ 1] > 0, 



when q[ 



it is easy to see that the third order doesn't vanish only 
- (2) - 2q[ 1) . In this case we obtain 



16 

x sin(2</> (1) 



(78) 



and 

& = ^(e^)\^sH2^-^) (79) 

x ~ft{q (1 \q (1) ) + ft{2q {1 \-q (1) ) + ft{-q {1) M 1) ) ■ 

Unlike the second order corrections, the Casimir energy 
depends on the phases <jyj> of the Fourier components of 
the corrugation. 

VII. CONCLUSIONS 



B. Nontrivial third order effects 

In the previous example, the third order correction to 
the Casimir energy vanished. We will now see that, when 
the corrugated surface can be described by a combination 
of sinusoidal functions of different wavelengths, a nonvan- 
ishing correction may arise to this order, when the wave- 
lengths satisfy a certain condition. Moreover, unlike the 
second order correction, the 0(?7 3 ) contribution depends 
on the relative phase of the sinusoidal components of the 
corrugation. 

Although one could consider a rather general corruga- 
tion, for the sake of definiteness we shall deal with the 
particular case 

lj(x||) = e (1) sin^xi + </> (1) ) + e (2) sin(gf >x x + (2) ) . 

. (74) 

The second order correction to the vacuum energy is eas- 



In this paper we have extended the perturbative re- 
sults for the Casimir energy between slightly deformed 
mirrors up to the fourth order in the amplitude of the 
deformations. For simplicity we considered the case of a 
deformed mirror in front of a plane mirror, for a scalar 
field satisfying Dirichlet boundary conditions. The third 
and fourth order results are important to improve the 
accuracy of the (almost) analytic perturbative calcula- 
tions, to have an explicit way to evaluate the validity of 
the rather simpler second order results, and to provide a 
benchmark for complex numerical calculations. The re- 
sults could be generalized, for example, to the case of the 
electromagnetic field in a geometry with two deformed 
mirrors. 

As a side point of the perturbative calculations, we 
have also discussed the relation between different approx- 
imations usually considered to compute Casimir forces. 
For the geometry considered in this paper, the Casimir 
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energy is a functional of the shape of the surface ip = 
a + rj. This functional can be expanded in powers of rj/a, 
as we did here, assuming small deviations from a flat sur- 
face. When the additional assumption that the function 
77 is slowly varying is reliable, a resummation of the per- 
turbative series is possible, and approximating the form 
factors by their zero-momentum values, the final result 
coincides with the PFA. This is a non perturbative re- 
sult, valid for arbitrary amplitudes as long as the surface 
is gently curved. The PFA can be improved by expand- 
ing the form factors around fc|| = 0. If this expansion 
only involves even powers of the momentum, the higher 
order corrections can be written in terms of derivatives 
of the shape of the surface. This is the derivative ex- 
pansion for the Casimir energy we proposed in previous 
papers [6 . On the other hand, when the form factors 
contain nonanalytic terms, the corrections to the PFA 
will be nonlocal, as it indeed happens for scalar field sat- 
isfying Neumann boundary conditions for nonzero tem- 
perature [7]. We have also obtained explicit expressions 
for the Casimir energy under the opposite assumption, 
i.e. strongly varying surfaces. In this case, the final re- 
sult could be written in terms of nonlocal operators, as 



Finally, we have presented some explicit examples. On 
the one hand, we have evaluated numerically the results, 
up to the fourth order, for the Casimir energy for the case 
of a sinusoidal corrugation. The numerical results repro- 
duce the analytic results expected in the limits of low 
and high wavelengths. The fourth order term is particu- 
larly relevant when evaluating the Casimir energy in the 
limit of short wavelengths, where it becomes dominant. 
The results presented in the paper allowed us to assess 
the validity of the second order calculations obtained in 
previous works, without the necessity of a full numerical 
evaluation of the Casimir energy. 



On the other hand, we have shown the case in which 
the corrugation consists of a combination of two sinu- 
soidal functions of different wavelengths, one being twice 
the other. The distinctive characteristic of the result for 
this case is that, when relation between the two wave- 
lengths holds, the Casimir energy bears a dependence on 
the relative phase of the sinusoidal functions appearing 
in the combination. This 'interference effect was absent 
in the second order term, namely, the relative phases of 
the components are irrelevant. 



APPENDIX 

In this Appendix we prove that the first terms in the 
perturbative expansion reproduces the PFA approxima- 
tion when the form factors are evaluated at zero momen- 



tum. The first term in the series of Eq. ( 63 1 is given in 



Eq.(39). Indeed, comparing both equations we obtain 



ftW(O) =a 4 



d 3 p\\ 



7T 

480 



(80) 



that reproduces the linear term in r\ in Eq. ( 65 1 



In order to evaluate the term quadratic in 77, we com- 
pare Eqs.([63} and (Eo). We get 



h (2) (0,0) = -a 5 



d?V\\ 



B(p||)(l+fl(p||))|p||| 



giving the quadratic term in Eq. ( 63 ) 



7T 

240 ' 
(81) 



The zero momentum contribution to the third order 



can be read from Eqs.(|43|), (|44|), and (|46|) 

/i (3) (0,0,0) =a 6 



(2ny 



hi 



B 

0~ 



B + B (4B 3 
2 3 V 



65 + 3) 



(82) 



where B = B(p»). The three terms in the integrals are 
the contributions coming from r( 3 J\ j = 1,2,3, respec- 
tively. The evaluation of the integral gives h^ 3 ' (0, 0, 0) = 
7r 2 /144, reproducing the third term in Eq.(65). 

The evaluation of the fourth order is more cumbersome 
but straightforward. The form factor (0, 0, 0, 0) can 
be written as a linear combination of f^ 4 ' l \ i = 1, ...5, all 
of them evaluated at zero momenta. Adding all contri- 
butions we obtain 



h (4) (0,0,0,0) = -a 1 



-B 



(2tt) 3 
7 



hi I 



B + AB 



2B 4 



(83) 



Evaluating the integral we obtain (0, 0, 0, 0) = 
— 7r 2 /96, in agreement with the fourth order term in 
Eq.pb. 
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